Evidence of transient interactions between it — > tt* optical excitations and image 

potential states in graphite 
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Here we report the experimental evidence of the interactions between the excitations of the it — > tt* 
optical transition and the image potential states (IPS) of highly oriented pyrolitic graphite (HOPG). 
By using non-linear angle resolved photoelectron spectroscopy (NL-ARPES) we show that the IPS 
photoemission intensity, the effective mass, and the linewidth exhibit a strong variation when the 
photon energy is tuned across the tt — > n* saddle points in the 3.1 - 4.5 eV photon energy range. 
A model based on the self-energy formalism is proposed to correlate the effective mass and the 
linewidth variations to transient many body effects, when a high carriers density (in the 10 20 cm -3 
range) is created by the absorption of a coherent light pulse. This finding brings a clear evidence 
of a high IPS-bulk coupling in graphite and opens the way for exploiting the IPS as a sensitive, 
nonperturbing probe for the many-body dynamics in materials. 

PACS numbers: 81.05.uf - 73.20.-r - 79.60.-i - 71.20.Gj 
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The concept of quasiparticlc in solids has been very 
useful for describing many unconventional physical prop- 
erties, such as the deviation from the linear response in 
many body systems, by taking into account the residual 
interactions between the quasiparticles, i.e. the part of 
the Coulomb forces not accounted for in the formation of 
quasiparticle 

For example, in these last decades it has been demon- 
strated that the optical properties of semiconductors near 
the electronic bandgap are dominated by excitons and 
their residual interactions, producing many-body effects, 
such as band-gap renormalization (BGR) and non- 
linear optical mechanisms. In particular, a deep insight 
has been gained on understanding the BGR in semi- 
conductors by means of time-resolved non-linear opti- 
cal spectroscopy in the fs time-domain [3(. For these 
reasons and also because their linear properties are well 
understood and excitons can be easily created and con- 
trolled through band-gap optical transitions, semicon- 
ductors have been considered ideal materials for such 
studies. Instead, in semi-metallic and gap-less systems 
the response to photoexcitation processes is still unclear, 
and the studies of the dynamics and interactions of quasi- 
particles have been only marginal. 

In a previous work Q , we have shown that for nonlin- 
ear photoemission processes in metals, the IPS effective 
mass (to*) can be influenced by the hot electron (hot-e~) 
gas generated during the excitation process, but the in- 
terplay between bulk many body excitations and surface 
states remains elusive. Nonetheless, understanding these 
processes is of paramount importance for condensed mat- 
ter physics. 

In this framework, graphite is an intriguing and inter- 
esting system, since the excitations in this material could 



have a strong coupling with the IPS, originating from the 
hybridization between the IPS and the interlayer band 
[5j , along with the anisotropy of the charge transport in 
the direction normal to the graphene planes. In addition, 
the presence of a van Hove singularity (vHs) at Hlj = 4.4 
eV in the JDOS, corresponding to the tt — > tt* transition 
between the bands saddle points, located at the M point 
of the Brillouin Zone (BZ) Q, allows the photoinjection 
of high carrier densities in the outermost graphene layers. 

Here we report on angle-resolved non- linear photoemis- 
sion spectroscopy (NL-ARPES) studies, using ultrashort 
coherent photon pulses, on HOPG. By tuning the photon 
energy in the 3.1-4.5 eV range we discover that the pho- 
toemission (PE) intensity of the n=l IPS, populated by 
the hot-e~ scattering process accompanying the tt — > tt* 
femtosecond excitations, has a maximum near ftw = 4.0 
eV. This energy results w 0.4 eV lower than the max- 
imum of the imaginary part of the dielectric function, 
e2(w), where the highest photoelectron yield is expected 
for indirectly populated IPS. 

Stimulated by this result and its possible link to tran- 
sient BGR effects, we measured also the IPS effective 
mass (to*) and linewidth (7tot), finding that both quanti- 
ties reach a maximum between huj = 4.0 eV and fko = 4.1 
eV. On the basis of the measured to* and 7tot versus the 
photon energy we craft a semi-quantitative self-energy 
model, allowing to obtain an estimation of the renormal- 
ization energy of « —0.4 eV. This value is fully consistent 
with a BGR mechanism that, by downward shifting the 
vHs energy, sets the IPS PE maximum in the hui w 4.0 eV 
range, in agreement with the experimental observations. 

The experimental setup is based on an amplified 
Ti:Sapphirc laser system and optical parametric ampli- 
fiers, producing p-polarized, ps 100 fs, near-UV pulses 
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that are focused on the HOPG sample kept in an ultra- 
high vacuum chamber (pressure < 2.0T0 -10 mbar). Pho- 
toelectrons are detected by a custom-built time of flight 
(ToF) electron spectrometer Q with an angular accep- 
tance of ±0.85° and an overall energy resolution 35 meV 
at an electron kinetic energy (Ek) of w 2 eV. The HOPG 
samples are cleaved ex- situ and annealed at 450° C until 
a good LEED pattern is obtained. The photoelectron 
spectra have been collected either at k|| =0 (normal 
emission - NE) and k|| + 1 0, being kn the component 
of the momentum parallel to the surface, by varying the 
angle between the sample normal and the ToF axis. In 
order to perform dispersion measurements with the max- 
imal axial symmetry possible, the plane of incidence of 
the incoming UV pulses has been chosen to lie, at NE, on 
the plane defined by the manipulator axis and the ToF 
axis, [see lower part of Fig. [lj-(a)]. All measurements 
have been performed at room temperature. 
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FIG. 1: (a) IPS population-PE process (red) and 7r band 
transition (blue) at fia;=3.14 eV (far from the vHs resonance, 
indicated by the empty blue ellipses). Solid (dashed) lines 
represents direct (indirect) photoinduced transitions. Gray 
wavy lines represent momentum-conserving scattering. Gray- 
filled, bell-like curves represent hot-e~ distributions; a sketch 
of the symmetric experimental geometry is shown. The inci- 
dent pulse is steered onto the sample at an azimuthal angle 
<j> = 45° with respect to the rotation plane (RP) and in-plane 
with the ToF and manipulator (MA) axis, (b) The same as 
for (a) at fiw=4.04 eV (near the vHs resonance, full blue el- 
lipses). The contribution to the IPS population from hot-e~ 
scattering from the M point is shown; the fc 2 = BZ is shown 
in the lower part of the figure; (k is the crystal momentum). 
Bulk graphite bandstructure is taken from Ref. @; Ef and 
Ey are the Fermi and vacuum energies, respectively. 



The upper parts of Fig. [JJ show normal emission NL- 
PE spectra of HOPG at hw = 3.14 eV (a) and hw = 



4.04 eV (b). The peaks located at Ek above 2 eV are 
identified with the n=l and n=2 IPS emission, instead 
the hump-shaped feature (tt*), located below 2 eV in Fig 
[JJ (b), is assigned to the linear PE from the tt* band, 
transiently populated via the tt — > tt* transition. These 
assignments are consistent with the s-polarized photon 
spectra (not shown here) where the IPS emissions are 
quenched and the tt* intensity is significantly reduced, 
according to the electric dipole selection rules for IPS 
and bulk bands with tt symmetry 

To properly describe these complicated non-linear pro- 
cesses, important insights can be achieved by measuring 
the PE multiphotonic order (MPO) at different photon 
energies, defined as the logarithmic slope of the IPS emis- 
sion peak intensity versus the pump fluence. We find a 
MPO=3 for hw = 3.14 eV and MPO=2 for hw = 4.04 
eV. These informations, along with the graphite band- 
structure [§], reported in the lower parts of Fig. [TJ shows 
that, at hw = 3.14 eV, the IPS population process is the 
result of a two-photon absorption mechanism, whereas 
at fiw=4.04 eV the IPS is populated by one-photon ab- 
sorption process. At both photon energies the IPS is in- 
directly populated trough electron scattering of the hot 
electron (hot-e - ) population created in the tt* band by 
the optical absorption (see Fig. [TJ. 

To unambiguously clarify the IPS population mecha- 
nism several NE PE spectra, a selection of which is shown 
in[5Ja), have been measured at a constant laser fluence F 
= 150 /iJ cm~ 2 by tuning the photon energy from hw — 
3.2 eV to hw = 4.2 eV. The MPO of the IPS photoemis- 
sion, versus photon energy, are shown in Fig. [2jb). As it 
is possible to note the MPO changes from 3 to 2 in going 
form 3.2 eV to 4.15 eV photon energy, being the transi- 
tion region ks 300 meV-wide and centered at hw — 3.9 eV. 
The linearity of the photoemission process from the IPS 
is shown in Fig. [2^c) where the IPS emission intensity is 
reported in a false colors plot in the (Ek, hw) plane. It is 
manifest that the relation between the Ek values of the 
IPS maxima and the photon energy is linear and it is well 
fitted by E K = hw + E B (black line), being E B the IPS 
binding energy referred to the vacuum level, Ey. Note- 
worthy, the IPS Eb value obtained from the data fit, i.e. 
Eb = —0.85 ±0.1 eV is in agreement with the only value 
found in literature (loj . These findings confirm that at 
hw < 3.8 eV an indirect two-photon population process 
plus a linear photoemission process (2+1=3 MPO) take 
place, whereas at hw > 4 eV the IPS photoemission spec- 
tra result from an indirect one-photon population plus a 
linear photoemission process (1+1=2 MPO). For further 
details, see Ref. I 111 . 

Being the population process of the IPS the result of 
scattering of the hot-e~ gas excited into the tt* bands, 
the IPS emission intensity should be proportional to the 
number of optically excited carriers in the tt* band when 
the photon energy is tuned through the tt — > tt* transi- 
tion, i.e. proportional to €2(w). From the NE PE spectra 
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FIG. 2: (a) NL-PE spectra of HOPG. A selection of spectra 
collected with photon energies between 3.65 eV and 4.15 eV is 
shown here, (b) n=l IPS MPO versus photon energy, showing 
a MPO=3 to MPO=2 transition in the 3.8-4.0 eV region, (c) 
The IPS region of the NL-PE spectra shown in (a) versus 
photon energy in a false-color representation. Black line is a 
linear fit of the IPS position versus photon energy, (d) IPS 
intensity (blue squares) versus photon energy. The solid line 
is a guide to the eye. The red dashed line represents the 
graphite £2 taken from ref. H3 . 



is possible to obtain the dependence of the IPS emission 
intensity versus the photon energy. Fig. |^d) shows that 
the IPS photoemission intensity depends on the photon 
energy, presenting a maximum at hoj ~ 4.0 eV. However, 
the £2(^)1 as taken from Ref. H2I has a maximum, re- 
sulting from the 7r — > n* optical transition between the 
two saddle points located at the M point of the Brillouin 
zone (BZ), i.e. at the vHs, located at fiui = 4.4 eV. From 
these data it is unambiguous that the IPS maximum is 
located « 0.4 eV below the £2(0;) maximum. The MPO 
changes, along with the graphite band structure, cannot 
explain this behavior, and, in particular, the sudden drop 
in the PE yield detected at Tilo > 4 eV. Moreover, no en- 
hancements of the linear photoemission, due to final state 
effects or other band structure mechanisms, are expected 
in the 4 eV range, from bandstructure analysis Q. 

A possible explanation for this discrepancy can be 
gained by considering a transient shrinking of the gap 
between the two saddle points originating from the high 
density of photoinduced carriers, which at an absorbed 
fluence of (1 - R)F w 100 cm~ 2 , are « 2 • 10 20 cm" 3 , 
i.e. two orders of magnitude greater than the room tem- 
perature carrier density at the Fermi level in graphite 
|13j |. Eq. EH has been used to estimate carrier densities, 
and a constant reflectivity R = 0.3 for the 3.2-4.5 eV 



range [14| has been considered. It is worth noting that 
these excitations are in the near-UV region and across 
the non-dispersive 7r bands saddle point, where the band- 
induced effective mass cannot be properly defined, and 
the contribution of the residual many-body interactions 
between electrons and holes should become important in 
"dressing" the quasiparticle. For these reasons, there are 
no theoretical predictions in the literature for such a case. 

Nonetheless, given the proposed role of the hot-e~ ex- 
citations to the IPS population processes, such a large 
BGR shift should imply a significant many-body contri- 
bution to the IPS properties, such as a deviation of the 
IPS effective mass and linewidth from the free-electron 
values. 

To verify this assumption, the IPS m* at photon en- 
ergies between fiw=3.14 eV and hu=4.50 eV have been 
measured from k| | dispersions, according to standard pro- 
cedures. The IPS dispersions are shown in Fig. [3fa), 
while the resulting m* versus photon energy is reported 
in Fig. H». In Fig. E^c) the IPS linewidth at k M = ver- 
sus photon energy, obtained from the NE spectra of Fig. 
|2fa), is shown. The m* versus the photon energy shows 
a maximum located at fko «4.1 eV, and corresponding to 
an increase of about 34% of the bare electron mass (m e ), 
as measured at hw «3.14 eV. 

Noteworthy, the effective mass rising is accompanied 
by a sudden increase of the IPS NE linewidth, occurring 
near 4 eV, from the reported value of 70 ~ 130 meV 
[po| to «180 meV. These experimental evidences strongly 
point towards a relevant coupling between the surface 
IPS and the bulk optical excitations in graphite. 

The IPS m* and linewidth dependence on photon en- 
ergy can be rationalized on the basis of a semi-empirical 
self-energy formalism, by considering the interaction be- 
tween the IPS and the photoinduced hot-e~ gas in the 
7T* band as the cause of the observed deviations from the 
free-electron case. Accordingly, a simple choice for the 
functional form of the real part of the IPS self-energy 
E w (k|| ) explaining the observed frequency-dependent dis- 
persion modifications is the following, 



ReE aJ (k||) = N 2 (uj) [a-fr 
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where 



N(lo) = (1 - R)F(hc)- 1 e 2 (uJ - duj) 



(1) 



(2) 



is the non-equilibrium photoexcited carrier density in the 
7T bands 15j; here cLo is the renormalization energy re- 
lated to the BGR shift, a and (3 are two parameters of 
the model: N 2 (ui)a represents the rigid shift of the dis- 
persion, whereas N 2 (w)f3 is the positive parameter rep- 
resenting the effective mass increase. From Eq. Q]is pos- 
sible to obtain the m* dependence on the photon energy 



m*(cj) =m c [l - (3N 2 (uj)\ 



(3) 
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FIG. 3: (a) IPS k|| dispersions at various photon energies 
(open circles). The parabolic fits (solid lines) are also shown. 
The free-electron parabola (dashed lines) is reported for com- 
parison. The hcu = 4.21 eV dispersion has been rigidly shifted 
(gray arrow) of about 0.1 eV upwards for clarity, (b) The IPS 
in* data extracted from dispersions are shown versus photon 
energy (open circles). The red solid line results from the fit 
of Eq. [3] (c) IPS PE linewidth versus photon energy (open 
triangles) from data reported in Fig. [2ja). The red solid line 
is the fit of Eq. H 



The contribution to the IPS line- width at ki|=0 from 
the scattering with hot-e~ can be obtained from the 
ImE^kii = 0) using Eq. [1] and the Kramers-Kronig 
(K-K) relations The total linewidth is given by, 



7totM =70 + aim- 
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where 7o=130 meV is the equilibrium IPS linewidth, and 
Tq = 0.07 eV is a Lorentzian broadening. Since ex- 
perimental data on dielectric constant are affected by a 
substantial experimental broadening, calculated values of 
e 2 (io) [l3| has been used. 

Fig. OJb) shows the comparison of the fitting of Eq. 
© (red curve) to the IPS m* (open circles), using the 
coupling constant (3 and the BRG shift dui as parame- 
ters. Instead, the fit of Eq. [4] to the linewidth data is 
reported in Fig. [3jc), using a and duo as parameters. The 
best dcu fit value results 0.4 ±0.1 eV from both m* and 
linewidth data. The fact that similar dco values are ob- 
tained from two independent measurements, i.e. m* and 
7tot, connected by the K-K relations and that this doj 
value is consistent with the energy shift between the IPS 
intensity and the 62(00) maximum as measured at equilib- 
rium, i.e. « 0.4 eV are considered as strong evidences of 
the photoinduced residual interactions between photoin- 



duced excitations and image potential states in graphite. 

Their strong interaction can be rationalized consider- 
ing the low interlayer coupling between graphene layers in 
graphite that prevents the diffusion of photoexcited car- 
riers in the bulk. Moreover, photoinduced excitations at 
the 7r-bands saddle point have already proven to strongly 
deviate from the Fermi-liquid behavior [lj| Il9j] ]. con- 
firming the peculiar character of the M point of the BZ 
in graphite. Studying the residual interactions depen- 
dence on the photodoping concentration would be im- 
portant in order to clarify the detailed properties of the 
many-body interactions. This would require performing 
fluence-dependent measurements over a substantial (i.e. 
at least a decade) fluence range. Unfortunately two con- 
curring causes such as space charge effects and low elec- 
tron statistics, limit the dynamic range of these experi- 
ments. However, it is possible to finely tune the absorbed 
fluence (and hence the it carrier concentration) by vary- 
ing the manipulator angle employing a non-symmetric 
experimental geometry, in which the plane of incidence of 
light is orthogonal to the manipulator axis; this induces 
angle-dependent variations of the absorbed fluence, re- 
sulting in a k| 1 -dependence in the measured IPS effective 
mass, in accordance with Eq. [3] (ioj . 

In conclusion, we have reported the experimental ev- 
idence of the interactions between the photoinduced ex- 
citations and the image potential states (IPS) on highly 
oriented pyrolitic graphite (EIOPG) when photoexcited 
by an intense (at F « 100 /iJ cm -2 ), ultrashort ( « 100 
fs) coherent light pulses, in the 3.1- 4.5 eV photon en- 
ergy range. The strong dependence on the photon energy 
of the IPS photoemission intensity, effective mass, and 
linewidth, observed by NL-ARPES, when frui approaches 
the 7r — > 7T* saddle points absorption has been attributed 
to many-body effects (residual interactions) originating 
from the high density of excitations. Finally, these find- 
ings bring a clear evidence of a high IPS-bulk coupling 
in graphite and opens the way for exploiting the IPS as 
a sensitive, nonperturbing probe for the many-body dy- 
namics in materials. 
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